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Abstract. We present evidence for variations in the fine-structure constant from 
Keck/HIRES spectra of 143 quasar absorption systems over the redshift range 0.2 < 
Zshs < 4.2. This includes 15 new systems, mostly at high-2 (z a b s > 1.8). Our most 
robust estimate is a weighted mean Aa/a = (-0.57 ±0.11) x 10 -5 . We respond to 
recent criticisms of the many-multiplet method used to extract these constraints. The 
most important potential systematic error at low-z is the possibility of very different 
Mg heavy isotope abundances in the absorption clouds and laboratory: higher abun- 
dances of 25 ' 26 Mg in the absorbers may explain the low-z results. Approximately equal 
mixes of 24 Mg and 25 ' 26 Mg are required. Observations of Galactic stars generally show 
lower 25,26 Mg isotope fractions at the low metallicities typifying the absorbers. Higher 
values can be achieved with an enhanced population of intermediate mass stars at high 
redshift, a possibility at odds with observed absorption system element abundances. 
At present, all observational evidence is consistent with the varying-a results. 

Another promising method to search for variation of fundamental constants involves 
comparing different atomic clocks. Here we calculate the dependence of nuclear mag- 
netic moments on quark masses and obtain limits on the variation of a and m q //lQCD 
from recent atomic clock experiments with hyperfine transitions in H, Rb, Cs, Hg + 
and an optical transition in Hg + . 



1 Introduction 



The last decade has seen the idea of varying fundamental constants receive un- 
precedented attention. The historical and modern theoretical motivations for 
varying constants, as well as the current experimental constraints, are reviewed 
in [1], the many articles in [2] and this proceedings. Here we set experimental con- 
straints on variations in two fundamental quantities, the fine-structure constant 
(a = e 2 /he; Sect. 2) and the ratio of quark masses to the quantum chromody- 
namic (QCD) scale (TO q //iQCD; Sect. 3), from optical quasar absorption spectra 
and laboratory atomic clocks respectively. 
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2 Varying a from Quasar Absorption Lines 

To serve the broad readership of this proceedings, we first outline the salient 
features of optical quasar absorption spectroscopy. Sect. 2.2 discusses in some 
detail the recently introduced many-multiplet (MM) method of constraining 
varying a from optical quasar absorption spectra. All optical data studied here 
were obtained at the Keck I 10-m telescope on Mauna Kea, Hawaii, with the 
High Resolution Echelle Spectrograph [3]. When applied to these data, this is 
the only method to date to yield internally robust evidence for a varying-a. We 
summarize our spectral analysis techniques and present this evidence in Sect. 2.3. 
After responding to some recent criticisms of the MM method in Sect. 2.4, we 
discuss the most important potential systematic error for the MM results - 
cosmological isotopic abundance variations - in Sect. 2.5. 

2.1 Quasar Absorption Lines 

The optical spectra of quasars are rich with absorption lines arising from gas 
clouds along the line of sight to Earth. For those unfamiliar with the anatomy 
of a quasar spectrum, Fig. 1 provides a tutorial. Of particular interest are those 
absorption clouds with sufficiently high hydrogen column density [number of 
absorbing atoms per area along the line of sight, iV(Hi)] to make metal absorp- 
tion lines detectable. These are classified as Lyman-limit and damped Lyman-a 
systems (LLSs and DLAs): LLSs have TV (Hi) > 2 x 10 17 cm _1 and DLAs have 
N(H i) > 2 x 10 20 cm _1 . It is the pattern of metal absorption lines - the relative 
separation between the different transitions - which carries information about 
the value of a in the clouds. The upper and lower panels of Fig. 1 detail some 
of the metal lines, those in the upper panel being of particular interest for the 
many-multiplet method in Sec. 2.2. Note the 'velocity structure' of the absorp- 
tion cloud: each transition comprises many 'velocity components', each of which 
probably corresponds to a separate absorption cloud, all components probably 
being associated with a single high redshift galaxy or dark matter halo. 

2.2 The Many-multiplet (MM) Method 

Initial attempts at constraining a- variability with quasar absorption spectra [5- 
10] used the alkali doublet (AD) method: for small variations in a, the relative 
wavelength separation between the transitions of an AD is proportional to a. 
While the AD method is simple, it is relatively insensitive to a-variations. The 
s ground state is most sensitive to changes in a (i.e. it has the largest relativis- 
tic corrections) but is common to both transitions of the AD. A more sensitive 
method is to compare transitions from different multiplets and/or atoms, al- 
lowing the ground states to constrain a, i.e. the many-multiplet (MM) method 
introduced in [11,12]. We summarize the advantages of the MM method in [13]. 

We first illustrate the MM method with a semi-empirical equation for the 
rclativistic correction, A, for a transition from the ground state with total angular 
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Fig. 1. Keck/HIRES spectrum of quasar GB 1759+7539 [4]. The full spectrum (middle 
panel) shows several broad emission lines (Ly-a, Ly-/3, Nv, Civ, Siiv) intrinsic to the 
quasar in the observer's rest-frame (i.e. vacuum heliocentric; bottom) and the quasar 
rest- frame (top). The dense 'Lyman-a forest' blue- wards of the Ly-a emission line is 
caused by cosmologically distributed low column-density hydrogen clouds along the 
line of sight to the quasar. The damped Lyman-a system (DLA) at z a b s = 2.625 and 
the Lyman-limit system (LSS) at z a b s = 2.910 give rise to heavy-element absorption 
lines red- wards of the Ly-a emission line (away from the confusing Ly-a forest). Some 
Zabs = 2.625 transitions are detailed in panels A-E & G-I. Even though the transitions 
in the top panels have very different line-strengths, the velocity structures clearly follow 
each other closely. Detection of many such transitions facilitates determination of the 
velocity structure and allows easy detection of random blends. For example, the blue 
portion of the Alii A1670 profile is blended with C iv A1548 in the z abs = 2.910 LLS 
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where Z is nuclear charge and many-body effects are described by C ~ 0.6. To 
obtain strong constraints on a- variability one can (a) compare transitions of light 
(Z ~ 10) atoms/ions with those of heavy (Z ~ 30) ones and/or (b) compare s-p 
and d-p transitions of heavy elements. For the latter, the relativistic corrections 
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will be of opposite sign which further increases sensitivity to a-variation and 
strengthens the MM method against systematic errors in the quasar spectra. 

In practice, we express the rest-frequency, u> z , for any transition observed in 
the quasar spectra at a redshift z, as 



where a z is a in the absorption cloud. For most metal transitions observed in 
quasar absorption spectra, the laboratory wavenumber, loq, is measured with low 
precision compared with that achievable from the quasar spectra (!) since, in the 
laboratory, the transitions fall in the UV. For example, despite a recent order 
of magnitude precision gain [14], the Civ A1548 and 1550 wavenumbers carry 
formal errors > 0.04 cm -1 . Compare this with the precision of «0.02cm _1 avail- 
able from absorption lines in a high resolution quasar spectrum (see Sect. 2.4). 
Dedicated laboratory measurements [15,16,14] of loq for many transitions now 
reach an accuracy of < 0.004 cm -1 allowing a precision of Aa/a ~ 10~ 7 to be 
achieved. Updated values of u>o are given in table 2 of [17]. 

The q coefficient of each transition contains all the relativistic corrections and 
measures the sensitivity of the transition frequency to changes in a. These have 
been calculated in [11,18-20] using many-body techniques. The accuracy of these 
calculations is given by how well various observable quantities (e.g. spectrum, 
g-factors etc.) of the ion in question are reproduced. To account for all dominant 
relativistic effects, the Dirac-Hartree-Fock approximation is used as a starting 
point. The accuracy is improved using many-body perturbation theory and/or 
the configuration-interaction method. For most transition combinations used in 
the MM method, the accuracy of these calculations is better than 10%. Note 
that in the absence of systematic effects in the quasar spectra, the form of (2) 
ensures one cannot infer a non-zero Aa/a due to errors in the q coefficients. The 
q coefficients used in this paper are compiled in table 2 of [17]. 

Figure 2 shows the distribution of q coefficients in (rest) wavelength space. 
Our sample conveniently divides into low- and high-z subsamples with very 
different properties (throughout this work we define z < 1.8 as low- z and z > 1.8 
as high-z). Note the simple arrangement for the low- z Mg/Feii systems: the Mg 
transitions are used as anchors against which the large, positive shifts in the Fe 11 
transitions can be measured. Compare this with the complex arrangement for 
the high- z systems: low-order distortions in the wavelength scale of the quasar 
spectra will have a varied and complex effect on Aa/a depending on which 
transitions are fitted in a given absorption system. This complexity at high-z 
may yield more robust estimates of Aa/a. The right panel quantifies this using 
simulations of the original 128 absorption system sample (see next section) which 
have been artificially compressed (see [17] for details). Even though the systems 
at high-z each respond differently to the compression of the wavelength scale, 
the binned plot reveals the average response is opposite to that in the low-z 
Mg/Fen systems. This is an important strength of the MM method: the low- 
and high-z samples respond differently to simple systematic errors due to their 
different arrangement of q coefficients in wavelength space. 
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Fig. 2. (Left) Distribution of q coefficients for transitions used in the MM method. 
For low-z Mg/Fen systems, a compression of the spectrum can mimic Aa/a < 0. 
However, the complex arrangement at high-z indicates resistance to such systematics. 
We define several 'g-types' by the horizontal bands and labels shown. (Right) Binned 
measurements of Aa/a from 20 simulations of 128 absorption systems (solid) and the 
same, real quasar absorption systems (dotted). For the top and middle panels we input 
the indicated values of Aa/a. The values and errors are recovered reliably. A wavelength 
compression is introduced for the bottom panel to reproduce the low-z quasar results. 
At high-z, the variety of q coefficients causes the expected large scatter but the average 
effect on Aa/a is opposite to that in the low-z systems. Simple distortions of the quasar 
spectra cannot explain the results 



2.3 Spectral Analysis and Updated Results 
Sample Definition 

Our data set comprises three samples of Keck/HIRES spectra, each observed in- 
dependently by different groups. The first sample [21] contains 27 low-z Mg/Fell 
systems. The second sample [22] contains transitions from a wide variety of ionic 
species (though mostly singly ionized; Al II, Al III, Si II, Cr II, Fe II, Ni II and Zn II) 
in 19 high-z DLAs and 3 low-z Mg/Fe n systems. An additional high-z DLA is 
from [4]. The third sample was graciously provided by W. L. W. Sargent and 
collaborators and comprises 78 absorption systems over a wide redshift range. 
Together, these samples comprise 128 absorption systems and form the total 
sample presented in [17]. To illustrate many points in the following sections, we 
provide example spectra of a low-z Mg/Fe n system and a high-z DLA in Fig. 3. 

In this work we update the second sample with 15 additional systems ob- 
served and reduced by two of the authors (JXP & AMW) and collaborators 
[23], containing a mix of low- and high-z systems. 

From the fully reduced spectra we select all systems which contain at least 2 
transitions of different g-type (defined in Fig. 2), thereby potentially providing 
a tight constraint on Aa/a. Only in cases where all selected transitions have 
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Fig. 3. Selected systems and transitions registered on a common (arbitrary) velocity 
scale. Data (histogram) are normalized by a fitted continuum. Our Voigt profile fit (solid 
curve) and residuals, normalized to the la errors (horizontal solid lines), are also shown. 
Tick-marks indicate individual velocity components. (Top) z a b s = f.2938 Mg/Fell 
system towards Q0636+680f . (Bottom) z abs = 2.9587 DLA towards Q1011+4315. 
Note that only optically thin components constrain Act /a strongly. 



very low signal-to-noise ratio (SN) do we not attempt a fit. Only in very high 
SN cases have we selected systems where only transitions of the same g-type are 
detected. Apart from the obvious issues of line-strength and possible random 
blends (Sect. 2.4), many well known instrumental limitations prevent us from 
detecting all MM transitions in every system. For example, the throughput of the 
telescope/spectrograph and detector sensitivity drops sharply below 4000 A and 
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above 6000 A. Typically, the spectrum is not recorded below 3500 A and above 
7000 A. Also, gaps in the wavelength coverage appear, particularly towards the 
red, since the spectrograph is an echelle-cross-disperser combination. Echelle 
orders cover ~ 60 A and inter-order gaps can be up to ~ 20 A. 

Profile Fitting 

For each system, the available transitions are fitted with multiple velocity com- 
ponent Voigt profiles. Each velocity component is described by three parameters: 
the absorption redshift z a b s , the Doppler broadening or b parameter and the col- 
umn density N. We reduce the number of free parameters by assuming either a 
completely turbulent or completely thermal broadening mechanism: correspond- 
ing components in all transitions have equal b or their 6s are related by the inverse 
square-root of the ion masses. To apply the MM method one must assume that 
corresponding velocity components in all fitted ions have the same redshift. This 
further reduces the number of free parameters. We discuss this assumption in 
Sect. 2.4. To each system a single extra parameter is added, Aa/a. This allows 
all velocity components to shift in concert according to their q coefficients. 

All free parameters are determined simultaneously using vpfit, a non-linear 
least-squares \ 2 reduction algorithm written specifically for analysis of quasar 
absorption spectra. The 1 a parameter uncertainties arc determined in the usual 
way from the diagonal terms of the final parameter covariance matrix. The as- 
sumption that off-diagonal terms are small (that parameters are not closely cor- 
related) is a good one for Aa/a: redshift and Aa/a are not correlated (Fig. 2). 
Monte Carlo simulations with 10000 realisations confirm the reliability of the 
parameter and error estimates [24] . It is important to realise that this numerical 
method ensures that constraints on Aa/a are derived in a natural way from op- 
tically thin lines and not from saturated ones. The derivatives of x 2 with respect 
to the saturated component rcdshifts are very small compared to the optically 
thin case and so only the optically thin lines strongly constrain Aa/a. If the two 
broadening mechanisms mentioned above result in significantly different Aa/a, 
the system is rejected. Otherwise, the broadening mechanism giving the lowest 
X 2 fit is selected. We also require that x 2 per degree of freedom, x 2 , is ~ 1. 

Results 

The distribution of Aa/a with redshift and look-back time as a fraction of the 
age of the Universe is shown in Fig. 4. We also provide basic statistics for the 
different samples and the total, raw sample as a whole in Table 1. Note that all 
three samples give consistent, significantly smaller values of a in the absorption 
clouds compared to the laboratory. Breaking the sample down into low- and high- 
z subsamples also yields consistent results despite the very different q coefficient 
combinations used (Fig. 2, left) and overall reaction to simple systematic errors 
(Fig. 2, right). We have conducted numerous internal consistency checks on these 
results, including direct tests of the wavelength calibration of the quasar spectra 
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Fig. 4. Aa/a and 1 a errors from the many-multiplet method for three Keck/HIRES 
samples: Churchill (hollow circles), Prochaska & Wolfe (triangles), Sargent (squares). 
Upper panel: unbinned individual values. Middle panel: binned results for each sample. 
Lower panel: binned over the whole sample. To calculate the fractional look-back time 
we use H = TOkms^Mpc- 1 , J? m = 0.3, Q A = 0.7, implying an age of 13.47 Gyr. 



Table 1. Statistics for different samples, xl is X 2 /(^s ys — 1) about the weighted mean 



Sample 


iVsys 




Aa/a (1Q- 5 ) 


xl 


Churchill 


27 


1.00 


-0.531 ±0.223 


1.109 


Prochaska & Wolfe 


38 


2.27 


-0.664 ±0.219 


2.024 


Sargent 


78 


1.76 


-0.620 ±0.129 


1.182 


Low-z (z < 1.8) 


77 


1.07 


-0.537 ±0.124 


1.064 


High-z (z > 1.8) 


66 


2.55 


-0.744 ±0.167 


1.739 


Raw total 


143 


1.75 


-0.611 ±0.100 


1.373 


FiduciaP 


143 


1.75 


0.573 ±0.113 


1.023 



Low-z sample + low-contrast sample ± high-contrast sample with increased errors. 



and the effect of removing individual transitions or entire ionic species from our 
fits. These are described in detail in [13,17]. 



Extra Scatter at High z 

Note that the scatter in the total low-z sample is consistent with that expected 
from the size of the error bars (i.e. \l ~ !)■ However, at high-z, Fig. 4 shows 
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significant extra scatter. This is reflected in the high \l values in Table 1. The 
weighted mean therefore exaggerates the true significance of Aa/a at high-z. 

We have identified the major source of this extra scatter at high z. Consider 
fitting two transitions with very different line-strengths (e.g. Alii A1670 and 
Nin A1709 in Fig. 3). Weak components near the high optical depth edges of 
the strong transition's profile are not necessary to obtain a good fit to the data. 
Even though the vpfit x 2 minimization ensures that constraints on Aa/a derive 
primarily from the optically thin velocity components, these weak components 
missing from the fit will cause small line shifts. The resulting shift in Aa/a is 
random from component to component and from system to system: the effect 
of missing components will be to increase the random scatter in the individual 
Aa/a values. This effect will be far larger in the high-z sample since only there 
do we fit transitions of such different line-strengths. 

We form a 'high-contrast' sample from the high-z (i.e. z a b s > 1-8) sample 
by selecting systems in which we fit both strong and weak lines, i.e. any of the 
Al II, Si II or Fe n transitions and any of the Cr II, Ni n or Zn n ones. To obtain a 
more robust estimate of the significance of Aa/a in this high-contrast sample, we 
have increased the individual 1 a errors until \t = 1 about the weighted mean. 
We achieve this by adding 1.75 x 10~ 5 in quadrature to the error bars of the 
27 relevant systems. Other procedures for estimating the significance at high-z 
are discussed in [17]. Fig. 5 identifies the high-contrast sample and presents the 
binned results with increased error bars. Table 1 presents the relevant statistics. 
The above procedure results in our most robust estimate from the 143 absorption 
systems over the redshift range 0.2 < z abs < 4.2: Aa/a = (-0.57 ±0.11) x 10~ 5 . 

2.4 Recent Criticisms of the MM Method 

Bckcnstcin [25] has pointed out that the form of the Dirac Hamiltonian for the 
hydrogen atom changes in a theory where a is considered a dynamical field. 
He warns that the resulting shifts in energy levels could be important for the 
MM method. Reference [26] have extended Bekenstein's model to many-electron 
atoms. They find that the energy shift within the Bekcnstein model is propor- 
tional to Aa/a defined in (2), up to corrections of order 1%. Thus, modifications 
to the energy shifts discussed by [25] for the hydrogen atom are not important 
for the MM method applied to many-electron metal ions. 

Bahcall et al. [27] criticized the MM method on many points, summarized 
in their table 3. Though none of these are real candidate explanations of our 
results, we address each criticism below to avoid future confusion: 

• Theory: Since one must calculate the q coefficients using sophisticated 
many-body techniques, [27] argue that the MM method is less reliable than, 
say, the AD method. The likely sources of error are discussed in detail by 
[20] and we give a flavour of them in Sect. 2.2. The q coefficients are known 
to high enough accuracy given our sample precision (discussed below). We 
again stress that if Aa/a is really zero, one cannot manufacture a non-zero 
value through errors in the q coefficients (2) if systematic errors in the quasar 
spectra are not important. 
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Fig. 5. The fiducial sample. Low- z sample (triangles), low-contrast sample (squares) 
and high-contrast sample with increased error bars (circles). The weighted mean, 
Act/a = (—0.57 ± 0.11) x 10~ 5 , is our most robust estimate for the HIRES spectra 

• Absolute or relative wavelengths?: Reference [27] argued that the MM 
method requires the measurement of absolute wavelengths in the quasar 
spectra 1 . This is incorrect. Since we simultaneously determine the redshifts of 
the absorption components and Aa/a (these are not degenerate parameters; 
Fig. 2) for each system, any velocity shift may be applied to the spectra and 
Aa/a will be unaffected. A velocity-space shift is a systematic error for the 
absorption redshifts, not for Aa/a. 

• Sample precision: It is trivial to estimate the expected precision available 
from our observational sample. We explained this calculation in [13] and re- 
peat it here. The Keck/HIRES pixels cover ~3kms _1 and so one reasonably 
expects to centroid barely resolved features (with SN ~ 30-50) to 0.3 km s -1 . 
One expects ~ 4 features per absorption system to be well-centroided in this 
way, providing a velocity precision of Av ~ 0.15 kms -1 or Auj w 0.02 cm -1 
for an ujq w 40000 cm -1 transition. The typical difference in q coefficient be- 
tween the Mgi/n and Fen lines is ^lOOOcm -1 and so, for a single Mg/Fell 
absorption system, (2) implies a precision of |Aa/a| ~ 1 x 10~ 5 . With ^50 
such systems, one expects a precision |Aa/a| ~ 0.14 x 10~ 5 (cf. Table 1). 




1 This suggestion is present only in preprint versions 1 & 2 of [27] (astro-ph/0301507). 
It is removed from later versions. Despite this, we address this point here to avoid 
further confusion in the literature. 
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As shown in Fig. 2 (and figs. A. 2 & A. 3 in [24]), simulations also provide a 
simple "reality-check" on the sample precision. 

• Line misidentification and blending: Reference [27] argue we may have 
misidentified many absorption features. In high resolution (R ~ 50000) spec- 
tra we largely resolve the velocity structure of absorption systems. Misiden- 
tifying transitions is highly improbable since, even by eye, the profiles of dif- 
ferent species follow each other to within |Aw| < lkms -1 . Confirming this, 
we obtain good fits to the absorption profiles with the number of free param- 
eters restricted by physical considerations (Sect. 2.3). Detecting blends from 
absorption at other redshifts is also greatly facilitated by high resolution. 
See [13,17] for thorough discussions of blending and Fig. 1 for an example. 
Even if we misidentified a small number of transitions in our sample and they 
miraculously mimicked the velocity structure of other detected transitions 
(thereby allowing a good fit), this would have a random, not systematic, 
effect on Aa/a. Indeed, compared with the AD method, the MM method 
is distinctly robust against misidcntifications and blends: many transitions 
constrain the velocity structure so identifying blends and misidentifications 
is all the more trivial. Also, any blends that are not identified have a smaller 
effect on Aa/a since many other transitions contribute to the constraints. 

• Velocity structure: The MM method assumes that corresponding veloc- 
ity components in different ions have the same redshift. Most transitions 
used are from ionic species with very similar ionization potentials and so ab- 
sorption from these species arise should arise co-spatially. Consider a Mg n 
velocity component Wueshifted with respect to the corresponding Fe n com- 
ponent by some kinematic effect in the absorption cloud. Clearly, this mim- 
ics Aa/a < for that component. However, kinematic effects would equally 
well redshift the Mgn components. Thus, the effect on Aa/a is random from 
component to component and absorption system to system. This argument 
is misunderstood in [27]: they feel it is unlikely that such random effects 
"average out to an accuracy of 0.2 km s" 1 over a velocity range of more 
than 10 2 kms~ 1 ", the latter quantity referring to the total velocity extent 
of a typical absorption system. Inspecting Fig. 4, it is clear that any extra 
scatter in Aa/a from kinematic effects derives only from the gas properties 
on velocity scales less than typical b parameters, i.e. < 5kms _1 . That we 
obtain excellent agreement between the velocity component redshifts in dif- 
ferent species to a precision |Aw| ~ 0.3 km s -1 illustrates this. If kinematic 
effects were important, they would be most prominent in the low-z values of 
Aa/a from the Mg/Fe n systems, appearing as an extra scatter beyond that 
expected from the 1 a errors. This is not observed. We discuss kinematic 
effects in more detail in [28] and [17]. 



5 Isotopic Abundance Variations 



we discussed in [29,30] and emphasised in [17,31], the main possible system- 
z error for the low-z results is that relative isotopic abundances may differ 
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Fig. 6. (Left) Isotopic structures for relevant Mg/Si transitions from measurements 
[33,34,15] or calculations [32]. Zero velocity corresponds to the structure's centre 
of gravity. (Right) Sensitivity of Aa/a to isotopic abundance variations. We alter 
the heavy isotope abundances proportionately: the heavy element fraction for Mg is 
r3 26 = const, x (0.10 + 0.11) where the numbers in parentheses are the terrestrial 
isotopic fractions of 25 Mg and 26 Mg. Much higher relative 25 ' 26 Mg abundances in the 
absorbers can explain the low-z results but the high- z results (containing the Sin 
systems shown) are insensitive to 29 > 30 Si abundances 



between the absorption clouds and terrestrial environment/laboratory, as sum- 
marized in Fig. 6. The isotopic structures used for Mg transitions are discussed 
in [17] whereas those for Si are from calculations similar to [32]. To our knowl- 
edge, isotopic structures for the transitions of Fen are not known. However, 
these should be far more 'compact' than those of, say, Mg since the normal mass 
shift will be > 5 times smaller. The results in Figs. 4 & 5 were obtained by 
fitting the quasar spectra with terrestrial isotopic abundance ratios. For those 
systems where Mg lines and no Si lines are fitted, these results correspond to 
the large square on Fig. 6, vice versa for the large triangle. The Mg n and Si n 
systems approximately correspond to the low- and high-z samples respectively. 
If we change the heavy isotope abundances we note the marked change in Aa/a 
for the Mg systems and the comparative insensitivity for the Si systems. This 
is expected given the distribution of q coefficients in Fig. 2 and the diversity of 
transitions fitted at high-z (see [17] for further discussion). 

In previous works we argued that the heavy isotope fraction for Mg in our 
absorbers is likely to be significantly less than the terrestrial value. This is based 
on observations of low metallicity 2 , Z, stellar environments in our Galaxy [35,36] 
and theoretical models of Galactic chemical evolution [37,38] where significantly 
sub-solar heavy isotope fractions are observed/expected at the low Zs of our 

2 Metallicity is the relative metal abundance with respect to that in the Solar System 
environment. For LSSs and DLAs, log 10 Z typically ranges from —2.5 to —0.5. 
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absorption clouds. However, some stars in [35,36] and those in some globular 
clusters [39,40] are found to have super-solar values. 

At low Z, the heavy Mg isotopes are thought to be primarily produced by 
intermediate mass (IM; ~2-8M ) stars in their asymptotic giant branch (AGB) 
phase [41-43]. Reference [38] included a contribution from IM AGB stars in their 
chemical evolution model, finding sub-solar heavy Mg isotope abundances at low 
Z, consistent with the above observations. Recently, [44] noted that enhanced 
heavy Mg isotope abundances might be produced at low Z if one assumes an 
IM-enhanced stellar initial mass function (IMF) at high z. However, such an 
IMF seems incompatible with current constraints. See [45] for general discussion 
of the observational constraints on the IMF. More recently, [46] find Ar, S and 
O abundances consistent with a normal IMF in a high-z DLA. See also [47]. 
For example, an IM-enhanced IMF could produce vast amounts of Fe via type 
la supernovae [48], in disagreement with Galactic and DLA abundance stud- 
ies. Moreover, AGB enrichment levels are constrained by the relative element 
abundances in the absorption clouds. For example, AGB stars produce large 
amounts of N relative to other enrichment processes (e.g. supernovae types la 
& II; e.g. [49]). However, the abundance of N relative to H is very low in DLAs, 
i.e. 10 _3 8 -10 -1 ' 5 solar [50,51]. Overall, these points are a barrier to ad-hoc IMF 
changes such as those suggested by [44]. Though enhanced heavy Mg isotope 
fractions are a possible explanation of the low-z varying-a results, we again 
conclude that they are an unlikely one. 



3 Varying a and m q /AQ C D from Atomic Clocks 
3.1 Introduction 

The hypothetical unification of all interactions implies that variation of the elec- 
tromagnetic interaction constant a should be accompanied by variation of masses 
and the strong interaction constant. Specific predictions require a model. For ex- 
ample, the grand unification model discussed in [52] (see also [53,54]) predicts 
that the quantum chromodynamic (QCD) scale, Aqcd, is modified as 

M 34 §^ (3) 

^QCD a 

The variation of quark and electron masses in this model is given by 

^ ~ 70^ , (4) 
m a 

resulting in an estimate for the variation of the dimcnsionless ratio 

5(mM QCD ) ^ 35 8a 
(m/AQCD) a 

The large coefficients in these expressions are generic for grand unification mod- 
els, in which modifications come from high energy scales: they appear because 
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the running strong coupling constant and Higgs constants (related to mass) run 
faster than a. This means that if these models are correct the variation of masses 
and strong interaction may be easier to detect than the variation of a. 

For the strong interaction there is generally no direct relation between the 
coupling constants and observable quantities, unlike the case for the electroweak 
forces. Since one can only measure variations in dimensionlcss quantities, we 
must extract from the measurements constraints on variation of m q //lQCD, a 
dimensionlcss ratio, where m q is the quark mass (with dependence on the nor- 
malization point removed). A number of limits on variation of m q //lQCD have 
been obtained recently from consideration of Big Bang Nucleosynthesis, quasar 
absorption spectra and the Oklo natural nuclear reactor which was active about 
1.8 billion years ago [9,55-64]. 

Below we consider the limits which follow from laboratory atomic clock com- 
parisons. Laboratory limits with a time base of several years are especially sen- 
sitive to oscillatory variation of fundamental constants. A number of relevant 
measurements have been performed already and many more have been started 
or planned. The increase in precision is very fast. 



3.2 Nuclear Magnetic Moments, a and m q /AQCD 

As pointed out by [65] , measurements of the ratio of hyperfine structure inter- 
vals in different atoms are sensitive to variation of nuclear magnetic moments. 
The first rough estimates of the dependence of nuclear magnetic moments on 
TO q //lQCD and limits on time variation of this ratio were obtained in [55] . Using 
H, Cs and Hg+ measurements [66], [55] limited the variation of m q //lQCD to 
about 5 x 10 _13 yr _1 . Below we calculate the dependence of nuclear magnetic 
moments on TOcj/zIqcd and obtain the limits from recent atomic clock experi- 
ments with hyperfine transitions in H, Rb, Cs, Hg + and optical transitions in 
Hg + . It is convenient to assume that the strong interaction scale ^Iqcd does not 
vary, so we will speak about variation of masses. 

The hyperfine structure constant can be presented in the following form, 



const, x 



[a 2 F iel (Za)] 



(6) 



The factor in the first bracket is an atomic unit of energy. The second, 'electro- 
magnetic', bracket determines the dependence on a. An approximate expression 
for the relativistic correction factor (Casimir factor) for the s-wave electron is 



7( 4 7 2 - l) ' 



(7) 



where 7 = — (Za) 2 and Z is the nuclear charge. Variation of a leads to the 
following variation of F Ic i [66]: 



rel 



Fel 



K— , 

a 



(8) 
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~ 7 2 (4 7 2 - 1) • 1 J 

More accurate numerical many-body calculations [18] of the dependence of the 
hypcrfinc structure on a have shown that the coefficient K is slightly larger 
than that given by this formula. For Cs (Z=55) K =0.83 (instead of 0.74), for 
Rb #=0.34 (instead of 0.29) and for Hg+ if =2.28 (instead of 2.18). 

The last bracket in (6) contains the dimensionless nuclear magnetic moment, 
H, in nuclear magnetons (the nuclear magnetic moment M — \i x eh/2M p c) and 
the electron-proton mass ratio, m e /M p . We may also include a small correction 
due to the finite nuclear size. However, its contribution is insignificant. 

Recent experiments measured time dependence of the ratios of hypcrfinc 
structure intervals of 199 Hg+ and H [66], 133 Cs and 87 Rb [67], and the ratio 
of the optical frequency in Hg + and the 133 Cs hyperfine frequency [68]. In the 
ratio of two hyperfine structure constants for different atoms, time dependence 
may appear from the ratio of the factors F re i (depending on a) and the ratio of 
nuclear magnetic moments (depending on ra q //lQCD)- Magnetic moments in a 
single-particle approximation (one unpaired nucleon) are: 



fx = [g s + (2j - l) gi ] /2 (10) 



for j = 1 + 1/2 and 



/i= 2(7TT) Hfe + (2j + 3)si] (ii) 

for j = I — 1/2. Here, the orbital g- factors are gi = 1 for valence protons and 
gi = for valence neutrons. Present values of the spin g- factors, g s , are g p = 
5.586 and g n — —3.826 for the proton and neutron. These depend on m q //lQCD- 
The light quark masses are only about 1% of the nucleon mass [m q = (m u + 
Wd)/2 « 5MeV]. The nucleon magnetic moment remains finite in the chiral 
limit of m u = m<j = 0. Therefore, one may think that the corrections to g s 
due to the finite quark masses are very small. However, there is a mechanism 
which enhances the quark mass contribution: 7r-meson loop corrections to the 
nucleon magnetic moments which are proportional to the 7r-meson mass ~ 
^/mq/lQCD- m 7r =140 MeV is not so small. 

According to [69], the dependence of the nucleon (/-factors on the w- meson 
mass can be approximated as 

■'<"•> = ird!W (12) 

where a = 1.37 GeV" 1 , b = 0.452 GeV~ 2 for the proton and a = 1.85 GeV" 1 , 
b = 0.271 GeV~ 2 for the neutron. This leads to the following estimate: 

^ = -0.174*^ = -0.087^, (13) 
g p m T m q 

^i^_0.213^ = -0.107^. (14) 

5n m v TO q 
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Equations (10,11,13,14) give variation of nuclear magnetic moments. For the 
hydrogen nucleus (proton) , 

^ = ^ = -0.087^. (15) 

For 199 Hg we have the valence neutron (no orbital contribution), and so 

*H = __i = -0.107^ . (16) 
M 9n m q 

For 133 Cs we have the valence proton with j=7/2 and 1=4, giving 

^ = 0.22^=0.11^. (17) 

For 87 Rb we have the valence proton with j=3/2 and 1=1, resulting in 

^1 = -0.128^ = -0.064^ . (18) 

Deviation of the single-particle nuclear magnetic moment values from the mea- 
sured values is about 30% and so we have attempted to refine them. If we neglect 
the spin-orbit interaction, the total spin of nucleons is conserved. The magnetic 
moment of the nucleus changes due to the spin-spin interaction because the va- 
lence proton transfers part of its spin, (s z ), to core neutrons (transfer of spin 
from the valence proton to core protons does not change the magnetic mo- 
ment). In this approximation, g s = (1 — b)g p + bg n for the valence proton (or 
g s = (1 — b)g n + bg p for the valence neutron). We can use the coefficient b as 
a fitting parameter to reproduce nuclear magnetic moments exactly. The signs 
of g p and g n are opposite, therefore a small mixing, b ~ 0.1, is enough to elimi- 
nate the deviation of the theoretical value from the experimental one. Note also 
that it follows from (13,14) that ?>g p /g p ~ 5g n /gn- This produces an additional 
suppression of the mixing's effect, indicating that the actual accuracy of the 
single-particle approximation for the effect of the spin g-factor variation may 
be as good as 10%. Note, however, that we neglected variation of the mixing 
parameter b. This is difficult to estimate. 



3.3 Results 

We can now estimate the sensitivity of the ratio of the hyperfine transition 
frequencies to variations in m q //lQCD- For 199 Hg and hydrogen we have 

S[A(Hg)A4(H)] 8a 8[m q /A QCD ] 
H(Hg)/A(H)] ^ a U - U KMqcd] " 1 j 



Therefore, the measurement of the ratio of Hg and hydrogen hyperfine fre- 
quencies is practically insensitive to the variation of light quark masses and 
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the strong interaction. Measurements [66] constrain variations in the parameter 
a = a[m q /A QCD ]-°- 01 : 



1 da 
a dt 



< 3.6 x lCT^yr- 1 . 



(20) 



Other ratios of hyperfine frequencies are more sensitive to m q //lQCD- For 133 Cs 
and 87 Rb we have 



8[A(Cs)/A(Rb)] _ 8a SKMqcd] 
[A(Cs)/A(Rb)] ' a KMqcd] ' 



(21) 



Therefore, measurements [67] constrain variations in the parameter X = 
a 0AS [m q /A QC „] - 17 : 



1 AX 

x~dT 



= (0.2 ±7) x lO -16 ^ -1 . 



(22) 



Note that if the relation (5) is correct, variation of X would be dominated by 
variation of m q /AQcr>'- (5) would give X oc a 7 . 
For 133 Cs and H we have 



8[A(G S )/A(R)} _ 8a 8[m q /4 QCD ] 
[A(Cs)/A(R)} ■ a + ' K/yl QCD ] 



(23) 



Therefore, measurements [70,71] constrain variations of the parameter Xh = 
a ' 83 [m q //l QCD ] - 2 : 

(24) 



1 dX H 



X H dt 



< 5.5 x lO^yr" 1 . 



If we assume the relation (5), we would have Xh oc a 8 . 

The optical clock transition energy, E(Hg), at A = 282 nm in the Hg+ ion, 
can be presented as 



E(Kg) = const, x 



m c e 



Frel(Za) . 



(25) 



Note that the atomic unit of energy (first bracket) is cancelled out and so we 
do not consider its variation. Numerical calculation of the relative variation of 
E(Rg) has given [18] 



5£(Hg) 



-3.2 



8a 



(26) 



£(Hg) 

Variation of the ratio of the Cs hyperfine splitting A(Cs) to this optical transition 
energy is equal to 



SL4(Cs)/£(Hg)] _ 6 Q 5a + 8K/Aqcd] + Q ^ SK/Aqcd] 



[^(Cs)/£;(Hg)] 



a [m e /A Q cu} 



[m q /A QCD ] 



(27) 



Here we have taken into account that the proton mass M p oc Aqqd- The factor 
6.0 before 8a appeared from a 2 F re i in the Cs hyperfine constant (i.e. 2+0.83) 
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and a-dependence of E(Hg) (i.e. 3.2). Therefore, the results of [68] give the limit 
on variation of the parameter U = a 6 [m e //lQCD][TOq//lQCD] ' 1 : 

< 7 x lCT^yr- 1 . (28) 

If we assume the relation (5), we would have U oc a 45 . Note that we present 
such limits on |(da/di)/a| as illustrations only since they are strongly model- 
dependent. 

4 Conclusions 

We have presented evidence for a varying a based on many-multiplet measure- 
ments in 143 Kcck/HIRES quasar absorption systems covering the redshift range 

0. 2 < z abs < 4.2: Aa/a = (-0.57 ± 0.11) x 10~ 5 . Three independent observa- 
tional samples give consistent results. Moreover, the low- and high-z samples are 
also consistent, which cannot be explained by simple systematic errors (Fig. 2). 
Our results therefore seem internally robust. The possibility that the isotopic 
abundances are very different in the absorption clouds and the laboratory is a 
potentially important systematic effect. A high heavy isotope fraction for Mg 
(r^ g 26 w 0.5) compared with the terrestrial value (-T^g 26 « 0.21) may explain 
the \ow-z results (Fig. 6). However, observations of low-metallicity stars and 
Galactic chemical evolution (GCE) models suggest sub-solar values of r^ 26 in 
the quasar absorption systems. GCE models with a stellar initial mass function 
greatly enhanced at intermediate masses may produce large quantities of heavy 
Mg isotopes via asymptotic giant branch stars. However, such models disagree 
with the observed element abundances in quasar absorption systems. The high- 
z results are insensitive to the isotopic fraction of 29 ' 30 Si. However, we stress 
the need for theoretical calculations and laboratory measurements of isotopic 
structures for other elements/transitions observed in quasar absorption systems. 

Aside from a varying a, no explanation of our results currently exists which 
is consistent with the available observational evidence. The results can best be 
refuted with detailed many-multiplet analyses of quasar absorption spectra from 
different telescopes/instruments now available (e.g. VLT/UVES, Subaru/HDS). 

We have calculated the dependence of nuclear magnetic moments on quark 
masses. This leads to limits on possible variations in m q //lQCD from recent lab- 
oratory atomic clock experiments involving hypcrfinc transitions in H, Rb, Cs, 
Hg+ and an optical transition in Hg + . These limits can be compared with limits 
on a-variation within the context of grand unification theories. Unfortunately, 
this comparison is strongly model-dependent. See, for example, [72]. 
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